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Abstract 

We have performed a fully non-perturbative calculation of the thermal properties of a system 
of spin 1/2 fermions in 3D in the unitary regime. We have determined the critical temperature 
for the superfluid-normal phase transition. The thermodynamic behavior of this system presents a 
number of unexpected features, and we conclude that spin 1/2 fermions in the BCS-BEC crossover 
should be classified as a new type of superfluid. 



I. INTRODUCTION 



The unitary regime is the situation in which the scattering length a is much larger than 
the average inter-particle distance: n|a| 3 ^> 1, where n is the number density^, 3|. At zero 
temperature, systems in this regime are widely believed to be superfluid, with a coherence 
length and an inter-particle distance of comparable magnitude. Such zero temperature 
problem has been considered by several authorsj^,^, ^|. It was shown experimentally in 2002 
that these systems are (meta)stable, and they have been extensively studied ever sincejf]]. 
From the theoretical point of view, the typical treatment is based on an idea put forward by 
Eagles, Leggett and others jl]. Their approach assumes a BCS-like form for the many-body 
wave function, which is then used for all values of a. The main problem with this treatment 
is that, close to the unitary regime, the fraction of non-condensed pairs becomes of order 
one a], and so a mean field description becomes questionable (even including fluctuations). 

To determine the thermal properties of fermions in the unitary regime, we have placed 
the system on a 3D cubic spatial lattice, with periodic boundary conditions. Since the 
system under consideration is dilute, the interaction that captures the physical situation is 
a zero-range two-body interaction V(ri — r 2 ) = — gSfoi — r 2 ), with a momentum cut-off hk c . 
Details of the method can be found in Ref. [9(. 



II. RESULTS AND CONCLUSIONS 

The results of our simulations are partially summarized in Fig. 1 (see Ref. |9| for further 
details). The thermodynamic quantities we computed present a number of features that 
can be easily identified. First, as T — > the energy tends to the T = results obtained 
by other groups . This confirms those results, as the algorithms they used suffer from 
a sign problem, which is not the case in the present approach. Second, there is a low 
temperature regime and a high temperature regime, separated by what we argue is the 
critical temperature for the onset of superfluidity, which we estimated to be T c = 0.23(2). 

For T < T c , the T-dependence of the energy can be accounted for by the two types of 
elementary excitation expected for this system: boson-like Bogoliubov-Anderson phonons 
and fermion-like gapped Bogoliubov quasiparticles. Their contributions can be estimated 
assuming that the system is a Fermi superfluid at T = 0, with a compressibility and pairing 
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gap as determined from results in Ref. j^y]. 

The chemical potential fi is essentially constant for T < T c , a fact reminiscent of the 
behavior of an ideal Bose gas in the condensed phase, even though we know the system is 
strongly interacting and superfluid. This unexpected result implies lack of fermionic degrees 
of freedom at those temperatures. 

Universality of the unitary regime together with fi(T) = const, for T <T C implies that 

£( r) = 4^g), { g)=f. + ((0\ n=\ (1) 

which is the temperature dependence of an ideal Bose condensed gas. According to our 
results n = 5/2 to about 10%. 

Above T c the system is expected to become normal. The energy behaves like the energy 
of an ideal Fermi gas, plotted in Fig. 1 with a vertical offset. This is surprising, because 
the estimated pair- bre ; kmg tenure b T* * 0,5. ( 8 ee Ref, Q Q , im pM» g t hat 
for T c < T < T* there should be a noticeable fraction of non-condensed pairs. Our results 




FIG. 1: The total energy E(T) is shown with open circles for a 8 3 -lattice and triangles for a 
6 3 -lattice. The chemical potential fi(T) shown with squares for the 8 3 -lattice. The Bogoliubov- 
Anderson phonon and fermion quasiparticle contributions E p h +qp {T) are shown as a dashed line. 
The solid line is Ep g {T) — 0.6efN(1 — where Ep g {T) is the energy of a free Fermi gas. Upper 
left inset: condensate fraction a(T) , with circles for 10 3 , squares for 8 3 and triangles the 6 3 lattices 
respectively. The solid curve is a(T) = a(0)[l - (T/T c ) 3 / 2 ]. 
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show no hint of their presence in this temperature interval. 

n 

The condensate fraction a(T), as defined in Ref. p, and evaluated at r = L/2 pair 
separation, is shown in the inset of Fig. 1. a(T) defines the off-diagonal long range order of 
U. t wo- body dens ity matrix 0. Tneten^e dependence O^T) is again consistent 
with T c = 0.23(2). Moreover, as T — > we recover the T = results in Ref. ^]. The 
T-dependence of a(T) resembles that of an ideal Bose gas, which comes as a surprise as 
well. 

At resonance, universality allows us to estimate the entropy as S(T) = (^E(T) — 
fi(T)N(T))/T, see Fig. 2. The knowledge of S(T) allows to establish a temperature scale 
at unitarity. Extending the suggestion of Ref. [l2j, from known T in the BCS limit, one 
determines the corresponding S(T BC s)- By adiabatically tuning the system to the unitary 
regime one can use S(Tbcs) = S(T unitar y) to determine T in the unitary regime. Especially 
at very low temperature, our results show relatively large errors in S due to both statistical 
and finite size effects. 

Summarizing, we performed a fully non-perturbative calculation of the thermal properties 
of a system of spin 1/2 fermions at unitarity. We determined the critical temperature for 
superfluidity to be T c = 0.23(2). The thermodynamic behavior of this system presents a 




T[ Ep ] 

FIG. 2: Entropy per particle as a function of temperature. Our data in circles for the 8 3 lattice 
and in squares for the 6 3 lattice; free Fermi gas in full line (corresponding to extreme BCS limit). 
Arrows indicate possible adiabatic cooling or heating processes. 
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number of unexpected features, suggesting that spin 1/2 fermions in the BCS-BEC crossover 
qualify as a new type of superfluid. 
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